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Abstract 

We investigate the decoherence properties of a central system composed of two spins 1/2 in 
contact with a spin bath. The dynamical regime of the bath ranges from a fully integrable inte- 
grable limit to complete chaoticity. We show that the dynamical regime of the bath determines 
the efficiency of the decoherence process. For perturbative regimes, the integrable limit provides 
stronger decoherence, while in the strong coupling regime the chaotic limit becomes more efficient. 
We also show that the decoherence time behaves in a similar way. On the contrary, the rate of 
decay of magnitudes like linear entropy or fidelity does not depend on the dynamical regime of the 
bath. We interpret the latter results as due to a comparable complexity of the Hamiltonian for 
both the integrable and the fully chaotic limits. 
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I. INTRODUCTION 



Real quantum systems always interact with their environment. This interaction entails 
that the system, initially in a pure state, becomes entangled with the environment and 
decays into an incoherent mixture of several states. This phenomenon, called decoherence, 
is an essential feature of quantum mechanical systems. 

From a fundamental point of view, decoherence provides a theoretical basis for the 
quantum-classical transition emerging as a possible explanation of the quantum ori- 
gin of the classical world. From a practical point of view, it is a major obstacle for building 
a quantum computer since it can produce the loss of the quantum character of the com- 
puter. Therefore, a complete characterization of the decoherence process and its relation 
with the physical properties of the system and the environment, such as the strength of the 
system-bath interaction, characteristic times of the bath, or the presence of quantum phase 
transitions or quantum chaos, is needed for both fundamental and practical purposes. 

Connections between decoherence and quantum chaos have been previously studied. How- 
ever, a universal theory has not yet been found. One line of argument establishes a link 
between the decoherence process and the Loschmidt echo j^, claiming that for a quantum 
system with a classically chaotic Hamiltonian the rate at which the environment degrades 
information on the initial state becomes independent of the system-environment coupling 
strength Another point of view, in some way contrary to the former, but generally 
accepted, states that a chaotic bath leads to faster and stronger decoherence than an inte- 
grable one . One significative manifestation of this phenomenon is the dependence 

of the decoherence time, i. e. the time for which the initial correlations in the central sys- 
tem are lost due to decoherence, with the system-bath coupling strength A. Some authors 
have found that for regular baths decoherence rate is proportional to A^, while chaotic or 
unstable ones display a considerable weaker dependence on A j^, Q]- A numerical study 
over a quantum walker with a complex coin has shown that, though a chaotic and a reg- 
ular environments may not be distinguishable in the short-time evolution, the chaotic one 
continues to be effective over exponentially longer time scales, whereas the regular bath sat- 
urates much sooner Sj]. A similar study on the Dicke model at weak coupling shows that the 
entanglement is smaller if the system is initially in a regular orbit than if it is in an irregular 
one However, exceptions for this general behavior are well known [l^. It is also argued 
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that, when the system-bath interaction becomes extremely small, so that the perturbation 



theory is applicable, the regular bath leads to a faster decoherence than the chaotic one 
Numerical studies of many-spin system show that a chaotic bath generates stronger and 
faster decoherence than an integrable one for strong enough coupling. However, the result 
is opposite in the perturbative regime [12|]. 

In this paper we study the connection between decoherence and quantum chaos in a 



12 1 in order to test if 



many-body spin systems. We follow the methodology proposed in 
the conclusions obtained there are applicable to a broader class of spin systems and, thus, 
can be postulated as generic. We use a Hamiltonian for the bath that depends on many 
arbitrary, real and independent parameters, and whose dynamical regime is independent 
of the specific values of these parameters. The integrable limit is defined as a random 
realization of the XYZ Gaudin magnet 13|, characterized by the existence of as many 



integrals of motion as quantum degrees of freedom. The transition to a chaotic regime is 
modelled by a single control parameter interpolating between the integrable Hamiltonian 
and a fully chaotic one. In both limits, and along the whole transition, the complexity of 
the Hamiltonian, understood as the number of different relevant terms, remains comparable, 
contrary to most of the previously studied systems, for which the regular limit is represented 
by a simplified Hamiltonian. For example, in [8[ the integrable limit is characterized by an 



12| it is reached when the Hamiltonian 



independent evolution of each spin of the bath, and in 
of the bath reduces to a site-dependent magnetic field, with a negligible interaction between 
different spins. 

In this work, we show that, for a wide class of spin Hamiltonians, the integrable limit 
generates decoherence more efficiently if the system-bath coupling strength is small, while 
the chaotic limit becomes more efficient when the coupling is larger. We also show that 
this conclusion can be extended to the decoherence time. Nevertheless, the transition from 
integrability to chaos in terms of magnitudes related to decoherence is not so smooth as it is 
on spectral statistics. Moreover, the rate of decay of the fidelity and the linear entropy does 
not depend on the dynamical regime of the bath, contrary to what it is usually claimed. 

The paper is organized as follows. In Sec. II, we describe the model and analyze its 
dynamical regime by means of spectral statistics. In Sec. HI, we study the efficiency of 
the decoherence process, illustrating the connection between the dynamical regime of the 
bath and some characteristic measures of decoherence, like the non-diagonal elements of the 
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system reduced density matrix and the linear entropy. In Sec. IV we study a quantitative 
characterization of the chaoticity of the bath in a perturbative regime using the linear entropy 
and the Loschmidt echo. We also study the relation between the onset of chaos and the 
decoherence time. Finally, in Sec. V we summarize our results. 



II. THE MODEL 



We will consider a central system composed by two interacting spins 1/2, Si and S2, and 
a bath composed of a large number of 1/2 spins 1^ [l^. The central system and the bath 
evolve with the following Hamiltonian: 

H = Hs + HsB + Hb, (1) 

where Hs is the self-Hamiltonian of the system, Hsb the interaction between the system and 
the bath, and Hb the Hamiltonian of the bath. For if 5 and Hsb we use the Hamiltonians 

Hs = JSi • S2, (2) 

and 

Hsb = Si ■ ^ akh- (3) 

k 

The interaction between the central system and the bath is carried out by a single spin of 
the system Si; the other spin, S2, is affected by the bath indirectly, through its interaction 
with Si, governed by Hs- These kind of models are useful to describe, for example, the 
destruction of Kondo effect by decoherence |l^. 

For the bath Hamiltonian Hb we use an XYZ model with long range interactions 

Hb = J2',H„ (4) 

j 

where 

N 

H, = J2 ^^^IPI + B,kI]Il + C,jpi (5) 
k=i^j 

and {ej} are free parameters. 

With this generic Hamiltonian we will describe a complete transition from integrability 
to chaos, depending on the properties of matrices A, B and C. The integrable limit is 
obtained when the Hamiltonians ([5]) fulfilled the conditions of the XYZ Gaudin integrable 



model. In this limit the N x N matrices A, B y C are defined in terms set of arbitrary 
parameters {zj}, according to the following identities: 

1 + K sn^{zj — Zk) 



A, 



k 



S\l[Zj - Zk) 
.2/ 



I — K sn [Za — Zk) 



C. 



jk 



sn[Zj - Zk) 
cn{zj - Zk) dn{zj - Zk) 



iin{Zj - Zk) 

where sn(u) = sn(u, /t) is the Jacobi elliptic function of modulus k, < k < 1, and cn(u) 
and dn(-u) are related to sn(-u) by d sn(u)/du = cn(-u) dn(M). The XYZ Gaudin model 



can be solved exactly by Bethe ansatz 13|]. There as many independent Hamiltonians 
as quantum degrees of freedom and, with the definition ([6]) they commute among them- 
selves, [Hi, Hj] Wi, j = 1, . . . , N for arbitrary values of the parameters {zj}. Therefore, they 



constitute a complete set of integrals of motion |l6[ |. 

The transition from integrability to the fully chaotic limit is performed by a single- 
parametric perturbation of the matrices defined above. If the amplitude of the perturbation 
is small, the resulting Hamiltonian is close to integrability; for increasing values of the 
parameter, the Hamiltonian approaches the fully chaotic limit. Such a perturbation can be 
achieved with the following identities 

A'j^ = (cos a) Ajk + (sin a) 

B'jk = (cos a) Bjk + (sin a) ii^^, (7) 
C'ji^ = (cos a) Cjk + (sin a) R^,^, 

where R^/^, and are random antisymmetric matrices, and < a < 7r/2. Therefore, 
for a > the algebraic structure of the integrable system is lost, in a similar way as the 
geometric structure of a classical integrable system is broken when it is perturbed (see 
for a complete discussion about this definition of quantum integrability and its connection 
with spectral statistics). Note, however, that Hb conserves its complexity, remaining a truly 
XYZ model along the whole transition; none of the X, Y and Z terms becomes negligible 
in the integrable limit. In consequence, this model allows to study the infiuence of the 
dynamical regime of the bath in decoherence process independently of its complexity. 
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A. Spectral statistics of the bath 



The concept of quantum chaos still lacks a clear definition. Usually, a quantum system is 
said to be regular or chaotic depending on the statistical properties of its spectrum. Using 
a semiclassical approximation, Berry and Tabor showed than the statistical properties of 
the spectrum of a generic quantum integrable system are well described by an uncorrelated 
Poisson distribution [l^. On the other hand, Bohigas et al. conjectured that the statistical 
properties of the spectrum of a generic quantum chaotic system coincide with those of 
Random Matrix Theory [2^. Therefore, the statistical properties of the spectrum of a 
quantum system are considered as a main signature of chaos in quantum mechanics (for a 



recent review see 



2l|) 



For a general quantum system, the level density p{E) can be separated into a smooth 



part p{E) and a fiuctuating part p{E). The former depends on the specific properties of the 
Hamiltonian, while the latter is universal depending only on the dynamical regime of the 



system 211]. Therefore, in order to determine whether a quantum system is regular or chaotic 
from the statistical properties of its spectrum, it is necessary to extract the fiuctuating part 
of the density. This is due by means of a procedure called unfolding, which maps every 



energy level Ei to a dimensionless magnitude C. 



Q = N{E,), (8) 
where N{E) is the accumulated level density 

N{E) = / dxp{x). (9) 



This map can be done analytically in a few simple systems, like quantum billiards or Random 
Matrix Ensembles, but in general it is a difficult task. In this paper, we have performed 



the unfolding by approximating N{E) with a set of Chebyshev polynomials by means of a 
least-squares fit. 

The most simple and widely used spectral statistic is the nearest-neighbor spacing 
distribution P{s), i.e. the probability distribution of the nearest-neighbor spacing se- 
quence Si = Ci+i ~ Ci- Foi' ^ regular quantum system the distribution follows a Poisson 
P{s) = exp(— s), while for a quantum chaotic system it follows a Wigner distribution 
P(s) = (7r/2) sexp (— 7rs^/4). Note that in both cases (s) = 1. 
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FIG. 1: P(s) distributions for regular a = (upper panel), and chaotic a = 7r/2 (lower panel) 
limits of Hb- In both cases there are = 13 spins, and the histogram is built by collecting 50 
different cases, characterized by different sets of parameters {ej}- The dashed line corresponds to 
the Poisson distribution, while the dotted line represents the Wigner distribution. 



In Fig. [T]we show the P{s) distribution for a set of 50 different realizations of the regular 
{a = 0) and fully chaotic {a = vr/2) limits of bath Hamiltonian Hb with = 13 spins. The 
random matrices R"*", and R^ are defined in terms of Gaussian random variables with 
zero mean, and variance equal to 

a(R^) = a(R^) = .(R^) = ^(A) + ^(B) + a(C) _ ^^^^ 

3 

Each realization is obtained by choosing an independent set of {ej} parameters, by means 
of Gaussian random variables with zero mean and a = 1. For all realizations the set of {zj} 
parameters is fixed to Zj = 3.71a/j7^5 order to cover the whole period of the Jacobi 
elliptic functions, and the modulus of the Jacobi elliptic functions is fixed to k = 0.5. As 
can be seen in Fig. [1] the regular limit clearly follows the Poisson distribution, while the 
chaotic limit is perfectly described by the Wigner distribution. 

In order to quantify the degree of chaoticity of the bath as function of the parameter a 
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FIG. 2: Parameter rj in function of a for = 9 (squares), = 11 (circles) and = 13 (triangles). 



it is useful to calculate the following quantity 

/so 
ds (P(s) - Pwigner{s)) 

'I = , (11) 

/ ds (-Ppoisson('S) -fWigner('S)) 
J 

where sq = 0.472913 determines the first intersection of Poisson and Wigner distributions. 
This parameter transits from rj = 1 to rj = when the system moves from integrability to 
chaos. Therefore, the curve //(a) shows how fast or slow is this transition. In Fig. [2] we 
show the value of as a function of a for three different sizes of the bath, = 9, = 11 
and A^ = 13. For the three cases the system very fast approaches to chaos for small values 
of a. The transition is increasingly faster for the larger bath sizes. 

The P{s) distribution describes short-range correlations, since it measures the fluctua- 
tions in distances between consecutive levels. To properly determine the chaoticity of a 
quantum system, it is also necessary to study the long-range spectral correlations. There 
are several statistics to measure this long-range correlations. The most commonly used are 
T^2{L) and A3(L) statistics j2l||. In this paper, we will use instead the 5„ statistic, defined 



from the unfolded energies as [22| 



= Cn+1 - Cl - ^- (12) 

This statistic measures the fluctuations of the unfolded energy levels {Q} from their average 
value. In particular, we are interested in its power spectrum 

N 



pS 



N 



"y^^ 5n ex'p(—27rink / N) 



n=l 



(13) 
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FIG. 3: Power spectrum of bn statistic for a = (upper left panel), a = vr/lO (upper right panel), 
a = 7r/5 (lower left panel), and a = 7r/2 (lower right panel). Squares correspond to = 9; circles 
correspond to = 11; and triangles correspond to = 13. The theoretical value for integrable 
system is plotted with a dashed line, and the theoretical value for chaotic systems, with a solid 
line. 



221. This 



which is proportional to for regular systems, and Xjk for chaotic systems 
statistic is simple to compute, and it is more sensitive to the dynamical regime of the system 
than the -P(s) statistic (see [23| for a detailed discussion of this point). 

We show in Fig. [3] the power spectrum of (5„ statistic for a = 0, a = tt/IO, a = 7r/5 and 
a = 7r/2 in a double logarithmic scale. The x axis is plotted in function of uj^ = 2TTk/N, 
which ranges from u = to u = tt independently of the size of the spectrum, and therefore 
allows to look for finite size effects quite easily. Both, the regular and the chaotic limit, closely 
follow the theoretical lines, except in the low frequency region. This region is spoiled by the 
unfolding procedure (see, for example, 2J] for a complete discussion about the misleading 



effects due to the unfolding). Note also than the smallest frequency available is ooi = 27t/N, 
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and thus larger systems cover a wider range of frequencies. For intermediate values of a, the 
result is similar to the obtained with the rj parameter, i.e. the transition to chaos is fast, 
and for a fixed value of a larger sizes are more chaotic. However, it is also seen that the 
transition for the power spectrum of 5„ statistic appears to be slower than the description 
given in terms of the parameter rj. In Fig. [21 the parameter rj identifies the system as almost 
chaotic for a = vr/lO and = 13, whereas in Fig. [3] (upper right panel) it is clearly seen 
that the statistic is still far from the chaotic limit (lower right panel). These differences 
can also be seen for a = vr/S, but they might be due to the spurious effects of the unfolding 
procedure. 

In conclusion, the bath Hamiltonian as defined in the previous subsection develops a 
complete transition from integrability to chaos, manifested in both short-range and long- 
range spectral statistics. Moreover, this transition is smooth and monotonous with the 
parameter a. The transition is faster for larger systems sizes and therefore a normalization 
a{N) is required for the results to be independent of the size of the bath. Hence, Hb is a 
good candidate to study the connection between decoherence and quantum chaos in finite 
size spin bathes. In the thermodynamic limit, the transition from integrability to chaos may 
be sharp for a > 0. 



B. Density of states of the bath 

Prior to the analysis of the connection between chaos and decoherence using the bath 
Hamiltonian defined above, it is important to check that the parameter a modifies the 
chaotic properties of the system without altering in a significant way the density of states of 
the bath. Changes in the density of states may have an important influence on decoherence 



processes 251] ■ In Fig. Hlwe show the density of states of Hb for three different values of 
a, corresponding to the regular limit, the fully chaotic limit and an intermediate case. For 
these calculations, the set of parameters {ej} has not been chosen randomly, but according 
to Cj = cos {V^j) ■ It can be seen that there are no important differences between the three 
cases under study. Similar results are obtained for other values of a (not shown) . Therefore, 
we can conclude that a determines the degree of chaoticity of the spin bath, without altering 
significantly the density of states of the spin bath. However, it is important to note that the 
density of states of the bath is quite sensitive to the parameters of the model Hamiltonian 
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FIG. 4: Density of states of with a = (squares), a = vr/lO (circles) and a = ir/l (triangles). 

Hb, specially to the set {zj}. The use of these kind of Hamiltonians to study the transition 
from integrability to chaos requires a careful definition of the parameters. 

III. DECOHERENCE AND CHAOS: LONG TIME EVOLUTION AND DECO- 
HERENCE EFFICIENCY 



A. General description of time evolution 

We will study the time evolution of the system from an initial product state: 

l^o) = l^o) ® \xo) ■ 

\iPq) is the initial state of the central system, for which we chose the singlet state 

(ITi)-liT)) 



(14) 



V2 



(15) 



Ixo) is the initial state of the bath, for which we chose a statistical superposition of all basis 
states with random coefficients. 

Depending on the size of the bath we can compute the time evolution either by means of 
a numerical diagonalization of the whole Hamiltonian H, or by a Chebyshev's polynomial 
expansion of the time evolution operator 6|]. In the former case we can treat up to = 11 
spins in the bath finding the exact evolution of the system. For larger systems we will resort 
to the approximation of the time evolution operator in a controlled way. In this section we 
will treat up to = 11 spins, which is large enough for our purposes. In the next section 
we will enlarge the bath up to = 15 spins using the Chebyshev expansion method. 
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To obtain a quantitative measure of decoherence, we calculate the reduced density matrix 
of the system 

pit) = TTs\^it)){^it)\, (16) 

where the subindex B indicates a trace over all degrees of freedom of the bath. In particu- 
lar, we will analyze the diagonal elements of the density matrix and the non-diagonal term 

(Til p(t) lit). At t = 0, (Til Pit) ITi) = (itl Pit) ITi) = 1/2, (TTI p(t) ITT) = (ii| pit) |ii) = 

0, and (Ti|p(^) lit) = (iT|p(^) ITi) = ^1/2- The time evolution shows how the entangle- 
ment between the system and the bath destroys the initial correlations of the system. 

Another useful quantity to measure the decoherence is the lineal entropy Q = Trp^. The 
initial state of the system is a pure state and the density matrix is idempotent (p^ = p), 
thus r2 = 1. The decoherence induced by the bath transforms state of system into mixed 
state with fi(t) < 1. Lower values of Q imply greater efficiency of the decoherence process. 

B. Long time evolution of system density matrix 

In Fig. [5] we show the time evolution for pi2 = (Til P lit); (Til P ITi) and (ii|p|ii), as 
a function of the parameter a and the system-bath coupling strength a, for very long times 
(note that the time axis is displayed in logarithmic scale). In all the cases we have chosen 
ttk = a V/c, and = 11 spins for the bath. We can see three main interesting features. 
First of all, for very long times the elements of the density matrix of the system relax to 
equilibrium states, called pointer states, which are relatively unaffected by the interaction 
with the environment and thus survive to the decoherence process. Second, as expected, the 
larger the value of the coupling strength constant the greater the efficiency of the decoherence 
process. And third, for small values of a, the regular system seems to give rise to a more 
efficient decoherence than the chaotic one, whereas for a ^ 1, this behavior is reversed. 

In order to better understand how chaoticity determines the efficency of the bath we 
calculate the pointer states elements of the system density-matrix of the previous calculation. 
The results are shown in Fig. [61 The elements of the system density-matrix in the pointer 
states, were obtained by averaging the results for t > 10^, i. e. Ppomt = {pit))tyio^- 
can see that in the three cases the regular limit gives rise to a more efficent decoherence for 
shorter values of the system-bath coupling strength, whereas this behavior is reversed for 
greater values of a. In the case of the non-diagonal pu elements, a ~ 1 is enough to totally 
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FIG. 5: Time evolution for different elements of the system density matrix: (Til P I it) (upper 
panel), (Ti|p|Ti) (middle panel), and (ii|p|ii) (lower panel), for regular q = (black line; blue 
online) and chaotic q = 7r/2 (grey line; red online) limits, as a function of the system-bath coupling 
strength a. 

destroy the initial correlations between If],) and Hf), characteristics of the initial state ipo. 
For the diagonal elements, however, a ~ 1 seems to produce the most efficent decoherence, 
since for a > 1 the values of the system density-matrix in the pointer state come back to 
the initial value. In this last two cases, the change in the relation between decoherence and 
chaos is clearly seen: for a < 1, the regular limit gives rise to stronger decoherence, while 
for a > 1 the chaotic limit becomes a more efficient. 
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FIG. 6: Matrix elements of pointer states of the system density matrix: (||| p|it) (upper panel), 
(Til/o|Ti) (middle panel), and (lower panel), for regular a = (squares) and chaotic 

a = 7r/2 (triangles) limits, and for an intermediate value a = vr/lO (circles), as a function of the 
system-bath coupling strength a. 



Similar results are obtained for the linear entropy Q. Fig. [7] shows the pointer state 
values of Q for a = 0, a = it /lO and a = tt/2. It is clearly seen that for a < 1, the regular 
bath produces a stronger decoherence, whereas for a > 1 this behavior is reversed. We can 
conclude that the dynamical regime of the bath determines the efficiency of the decoherence 
process. Moreover, the system-bath coupling strength determines whether integrability or 
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FIG. 7: Pointer states values of the system linear entropy, for regular a = (squares) and chaotic 
a = 7r/2 (triangles) limits, and for an intermediate value a = tt/IO (squares), as a function of the 
system-bath coupling strength a. 

chaos give rise to more efficient decoherence. 

IV. DECOHERENCE TIME AND SHORT TIME EVOLUTION OF SYSTEM 
DENSITY MATRIX IN PERTURBATIVE REGIME 

Having established the relation between the efficiency of decoherence and the dynamical 
regime of the bath for different values of the system-bath coupling strength, we now focus 
on the quantitative analysis of how quantum chaos affects the decoherence process. Our 
aim is to determine whether a smooth transition from integrability to chaos in the bath can 
be detected from quantities directly related to decoherence. For this purpose, we select a 
small value for the system-bath coupling strength, for which the integrable limit produces 
a stronger decoherence than the chaotic one. This choice allows us to follow the central 
system decoherence by its fidelity F{t), analogous to the Loschmidt echo, which measures 
the sensitivity of the system to external perturbations. The fidelity is defined as 

F{t)^Trs[p\t)p{t)], (17) 

where p'(t) is the system density matrix for an ideal evolution in which the system and the 
bath do not interact, i. e. Hsb — 0. The subindex S denotes a trace over the states of the 
central system. This quantity behaves in a similar way as the linear entropy Q. 

In what follows, we will consider a bath composed of = 15 spins. The evolution is 
approximated by means a Chebyshev expansion of the evolution operator. We have checked 
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FIG. 8: Fidelity F{t) for the central system perturbed by a bath composed by = 15 spins with 
Ofc = Y^^IO^^ yk, for different values of parameter a: solid black (blue online) corresponds to 
a = 0; dashed black (green online) corresponds to a = 7r/20; dotted grey (cyan online) corresponds 
to a = vr/lO; dotted black (black online) corresponds to a = it/ 5; dashed grey (magenta online) 
corresponds to a = Svr/lO; and solid grey (orange online) corresponds to a = 7r/2. 



that the size of the bath does not change qualitatively the decoherence process, but the 
fluctuations around the pointer states (see Fig [5]) are decreased. The system-bath coupling 
strength is set to = y^lO~^ V/c. This particular value was chosen because the coupling 

behaves as 6 = \J^^^k=i^k therefore the set of {0^} has to be scaled with the bath size. 

In Fig. [8] we show the fidelity F{t) for different values of parameter a. Several interesting 
facts emerge from the figure. First of all, the shape of all the curves is very similar. The 
main differences between them are related to the pointer states, i. e. the values reached 
after a long time evolution. We can also see that these pointer values of the fidelity do not 
increase monotonically with parameter a, as it is expected since the integrable limit gives 
rise to a stronger decoherence than the chaotic one for this value of the system-bath coupling 
strength. The numerical results show that when a is close to zero and, thus, the bath is 
close to integrability, the decoherence is more efficient than when a is close to 7r/2, that is, 
when the bath is close to the fully chaotic limit. However, this efficiency does not decrease 
monotonically for increasing a. In fact, the curves in the figure show that for a = vr/20, the 
fidelity is smaller than for a = 0, and for a = Svr/lO, it is larger than for a = 7t/2. Secondly, 
the transition from the values characterizing integrability to those corresponding to chaos is 
slower than the transition from integrability to chaos determined by the spectral statistics. 
For example, for a = vr/lO, the fidelity is close to the integrable limit a = 0, whereas the 
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0.6200 


0.8781 


127.1 


7r/20 


0.6118 


0.9146 


121.0 


vr/lO 


0.6306 


0.9036 


121.0 


7r/5 


0.6817 


0.8691 


134.9 


Svr/lO 


0.7088 


0.8402 


148.2 


7r/2 


0.6797 


0.8589 


135.0 



TABLE I: Parameters of eq. (jl8p for different values of a. 

spectral statistics are closer to the chaotic limit. In particular, the parameter rj indicates an 
almost chaotic behavior for a = it/ 10 (note that from the results shown in Fig. [2] we can 
conclude that an increase of the size of the bath accelerates this transition, and therefore we 
may expect that for = 15 and a ~ 7r/20 the bath is almost chaotic). The power spectrum 
of 6n statistic is also close to the theoretical value for the fully chaotic limit, but it reveals 
that such limit is not yet reached. 

Another important quantity related to the decoherence process is the decoherence time, 
that is, the characteristic time for the loss of coherence of the central system due to the 
coupling with the bath. One way to estimate this time is by means of the decay that the 
bath induces in F{t). We can fit the shape of F{t) to the following expression 

F{t) =p+(l-p) exp (-t^^i") , (18) 

where p is the pointer value for F{t), and (3 and are free parameters, the last one cor- 
responding to the decay time of the system. The results for this fit are shown in table [B 
We can see that the decoherence time is slightly larger for a chaotic bath, and that there 
is a clear correlation between the pointer value p and the decay time T^. Moreover, table 
[I also displays a surprising result: contrary to what is expected, F{t) decays roughly in an 
exponential way for all the values of a. Therefore, it seems that the dynamical regime of 
the bath does not affect the decay ratio of the Fidelity. 

The close connection between the fidelity F{t) and the linear entropy Q{t) can be seen 
in Fig. [9l The shape of the curves for linear entropy and fidelity are almost identical. The 
only appreciable difference between these two quantities is their pointer value. Moreover, 
the transition from integrability to chaos follows the same trend: it is non-monotonic with 
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FIG. 9: Linear entropy ^}{t) for the central system perturbed by a bath composed by = 15 spins 
with Ofc = ^J^W^'^ yk, for different values of parameter a: solid black (blue online) corresponds to 
a = 0; dashed black (green online) corresponds to a = 7r/20; dotted grey (cyan online) corresponds 
to a = vr/lO; dotted black (black online) corresponds to a = it/ 5; dashed grey (magenta online) 
corresponds to a = Svr/lO; and solid grey (orange online) corresponds to a = 7r/2. 



a 


P 


/3 


Ts 





0.4522 


0.8801 


112.5 


7r/20 


0.4459 


0.9078 


106.5 


vr/lO 


0.4655 


0.8994 


107.7 


7r/5 


0.5228 


0.8796 


125.0 


Svr/lO 


0.5515 


0.8599 


140.6 


tt/2 


0.5138 


0.8700 


130.2 



TABLE IL Parameters of eq. (I18p applied to lineal entropy ^}{t) for different values of a. 

a, and slower than the corresponding transition in the spectral statistics. 

A characteristic time of decoherence can be also defined from Q{t) fitting the numerical 
results to an expression similar to f|T8|) . We show the results in table[Tll They are very similar 
to those obtained with the fidelity: for values of a closer to the integrable limit, decoherence 
takes place faster than for values closer to the fully chaotic regime. There is also a strong 
correlation between the pointer value and the decoherence time. We also remark that the 
transition is not monotonous with a. 

From all these results we can conclude that the dynamical regime of the bath influences 
the pointers states of the central system in a non-strictly monotonic way. The modifications 
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of the pointer states are not smooth nor monotonic when the bath change from integrabihty 
to chaos. We have also shown that the decoherence time is influenced by the dynamical 
regime of the bath in a very similar way: from the decay of both Loschmidt echo and 
linear entropy, a trend from integrabihty to chaos can be identified in the decoherence time. 
However, it is also important to note that the rate of this decay is not affected by the 
dynamical regime of the bath, since it is roughly exponential along the whole transition. 
A possible explanation of this surprising feature is that the complexity of the Hamiltonian 
remains more or less the same along the whole transition from integrabihty to chaos, because 
it always consists on a XYZ model in which none of the terms is negligible. 



V. CONCLUSIONS 



We have studied the decoherence of a two-spin central system interacting with a bath 
whose dynamical regime can transit from integrabihty to chaos. Unlike previous studies the 
integrable regime is described by an XYZ Gaudin magnet with random parameters, with 
a complexity similar as the corresponding to the fully chaotic regime. We have calculated 
time evolution by numerical technics, and we have analyzed several quantities related to the 
reduced density matrix of the central system p{t). 

From our results, we conclude that at t — > oo the reduced density matrix of the system 
p has an observable dependence on the dynamical regime of the bath. For small values of 
the system-bath coupling strength, the asymptotic value of (ti| p lit) is larger in the regular 
limit than in the chaotic limit. However, this difference tends to decrease as the system- 
bath coupling strength is increased. For the diagonal elements (TilpiTi) and (ii|p|ii), 
the relation between the dynamical regime of the bath and the pointer states changes at 
a ~ 1. Below this value, the regular limit gives rise to a stronger decoherence; above it, the 
efficiency of the decoherence process is larger when the bath is chaotic. For both diagonal 
elements a ^ 1 gives rise to the larger decoherence. These results show that the onset 
of chaos affects the decoherence process of the central system in a non simple way, since 
its influence depends on how strong is the coupling with the bath. These conclusions are 
consistent with those obtained in 12 1. 

The connection between the properties of the reduced density matrix of the system and 
the dynamical regime of the bath allows to analyze the transition form integrabihty to chaos 
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in terms of decoherence. We have done so with a small value of the system-bath coupling- 
strength, that is, in the perturbative regime. Our results show that this transition is not so 
smooth and monotonous as described by the spectral statistics analysis, in spite of the fact 
that the regular and chaotic limits are clearly distinguished. Beyond this non-monotonic 
behavior, we have also shown that the transition from integrability to chaos in the central 
system reduced density matrix is slower that than the transition in the nearest neighbor 
spacing distribution by means of rj parameter, and it is a bit closer to the behavior of the 
power spectrum of (5„ statistic. Therefore, long-range correlations in spectral fluctuations 
seem to be involved in the bath efficiency to produce decoherence. 

We have also performed a similar analysis with the decoherence time, that is, the time 
at which the central system losses its original correlations. For a perturbative regime, the 
integrable limit produces stronger decoherence in shorter times as compared with the fully 
chaotic regime. Nevertheless, in contradiction with the results of reference 12|, the decay 
rate of the fidelity and the linear entropy does not depend on the dynamical regime of 
the bath. The main difference between both treatments is the modelling of the integrable 
regime; while our integrable limit contains all the complexity of the chaotic regime, reference 



12l | uses a simple integrable limit with large degeneracies. Therefore, we conclude that the 



decay of the fidelity and the linear entropy is related to the complexity of the bath and not 
to its dynamical regime. Further work is needed to clarify this result. 
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